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Abstract 

The effect of a random magnetic field in the convective zone of the Sun on resonant 
neutrino spin oscillations, i.e. transitions of the type v e L — > ^fl, is considered. The 
average survival probability and the expected experimental signals in the existing 
solar neutrino experiments are computed as a function of the level of the noise and 
magnitude of a constant magnetic field in the convective zone. From comparison 
with observed detection rates we conclude that the RSFP solutions to the SNP with 
negligible mixing angle are stable under the presence of low or moderate levels of 
noise. Detection rates, specially in the Homestake experiment, are however sensitive 
to large levels of noise. As a consequence, an upper limit on small scale magnetic 
fluctuations is obtained from the combined solar data: yj (B 2 ) < 140 — 200 kG for the 
scale Lq r~j 1000 km and transition moment /i = 10~ 11 /j,b- 
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1. A neutrino transition magnetic moment can account, as it has been shown in 
for both the observed deficiency of the solar neutrino flux and the possible time variations 
of the signal. The overall deficit is caused by the alteration or suppression of the neutrino 
energy spectrum. The time dependence may be caused by time variations of the magnetic 
field in the convective zone of the Sun. 

The two main problems for the RSFP scenario are the large size of neutrino magnetic 
moment which is needed and the scarce knowledge which is available on the structure of the 
magnetic field inside the Sun. To overcome the second problem, authors usually consider 
different "plausible" radial magnetic field profiles. The absence of a strong poloidal field 
around the Sun limit strongly the maximum value of the magnitude of a magnetic field 
in its inner parts. Therefore, a zero magnetic field is usually considered for the solar core. 
In the convective zone, fields of the order of a few kG up to a few tens of kG are usually 
considered possible. 

In this article we will address another aspect of the RSFP scenario: the nature of 
the magnetic field. Magnetic fields in nature are very often stochastic or at least strongly 
chaotic. This happens both, when considering the magnetic fields coming from the seeding 
of primordial fields or more mundane fields such as those present in the external shell 
of the Sun. Neutrino propagation in a medium with random magnetic fields has special 
properties, differentiated from the situation where only regular magnetic fields are present. 
The effective Hamiltonian governing the evolution of averaged quantities become non- 
hermitic and some sort of dissipation is introduced. Average conversion probabilities 
become aperiodic and more sensitive to vacuum parameters in some cases. It has been 
shown for example that limits on magnetic moments derived including random magnetic 
fields in Supernova are stronger than those derived using only regular fields (see ||]). On 
the other hand, for strong noise the situation is in some sense opposite, as we will see 
below in this work the survival probability approaches one half irrespective of the value 
of the vacuum parameters. 

The aim of this letter is to study the effect of the introduction of a random magnetic 
component in a simple RSFP scenario: one which disregards neutrino mixing and assumes 
that the direction of the transverse field does not change along the neutrino trajectory. 
It has been shown (Q) that the same model with a regular magnetic field profile can 
provide a rough solution to the SNP. The conclusion of ||] is that all the available solar 
neutrino data can be fitted for Am 2 ~2x 1CP 8 — 4 x 10 -9 eV 2 with no inner field and a 
convective-zone magnetic field varying between 2.5 — 5 x 10 4 G (fj, = 10 -11 /j-b)- 

Summarily, the main advantages of models including random magnetic fields are con- 
nected with: a) The no neccesity to know in detail specific magnetic profiles, b) The 
potentially stronger effect of the random magnetic field in comparison to regular mag- 
netic field: for the latter the spin precession is roughly proportional to the quantity 
expC J drB{r) where the integral is extended over the region where the magnetic field B 
is present. However for the case of random magnetic fields the precession is proportional 
to expC'J drB 2 (r). In addition the chaotic precession is irreversible and regeneration ef- 
fects of the original neutrino flux are absent or reduced, c) It is still important to assume 
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a large transition magnetic moment // for the neutrino, but the dependence on it is sub- 
stantially altered. While for regular oscillations the constant C ~ [a, for random magnetic 
fields C ~ S7 2 Lo/i 2 . The dependence on the magnetic moment [i can be traded off, at 
least partially, for the dependence on the parameters Q 2 ,Lq. These parameters describe 
the randomness and the coherence length of the magnetic field, they can be potentially 
very large, with an accepted range of variation of ~ 1 — 3 orders of magnitude, d) The 
overall influence of the random nature of the magnetic field is the flattening of the neutrino 
spectrum. The strong energy dependence on resonances is reduced, so this scenario is a 
natural way to avoid strong time variations. 

The structure of the article is as follows. First we will write the basic neutrino prop- 
agation equations in presence of matter and magnetic field. In the next section a brief 
description of the corresponding evolution equation for the averaged density matrix (the 
Redfield equation) is given. Finally the last section is dedicated to show the results of the 
numerical calculations and concluding remarks: The Homestake, SAGE-GALLEX (Ga- 
Ge) and (Super)-Kamiokande expected signal rates has been calculated and compared 
with available experimental observations. 

2. Taking into account theoretical and observational evidence of the solar surface and 
convection region magnetic fields it is expected (||, ||, [f| 0, ||]) that the neutrino, while 
traveling through the convection zone, will encounter a magnetic field of complicated un- 
known structure which will be the sum of large and small (fluctuating) scale components of 
sufficiently large amplitude and short correlation length. In this work, we will approximate 
such a field by the simplest possible expression: 

B{t) = B + B(t) (1) 

where Bq is a constant magnetic field over the full extent of the convective zone. B{i) is a 
fluctuating field with zero average acting over a relative small layer at the bottom of the 
convective zone (the unstable transition region between r ~ 0.65 — 0.75 Rq). It is assumed 
that B{t) is a stochastic field: a zero-average 5-correlated Gaussian process characterized 
by a two-point correlation function parametrized as follows: 

(B(t)B(t')) = (B 2 )L 5(t-t') = V B 2 L 5(t-t'). (2) 

The length Lq is a basically unknown parameter, values in the range 10 3 - 10 4 should 
be allowed. It has been shown in || that the 5-correlation function is a sufficiently good 
approximation to more realistic finite correlators even for relative large correlation lengths. 
The coefficient rj is the ratio between large and small scale fields: it is estimated ( pj], p. 
517; see also ||) by the expression 

1 v fo\ 
e Vo 

where Vo corresponds to the Alfven speed in the large scale field Bq, v is the turbulent 
fluid velocity and e is a, model dependent, undetermined number smaller than one. In the 
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turbulent dynamo theory the magnetic field is carried in the plasma; equating magnetic 
and kinetic energy, we obtain v ~ Vo, and so rj > 1. Some other estimations favor much 
higher values ?] >> 1 (p0|). 

In transverse magnetic fields, neutrinos with transition magnetic moments will experi- 
ence spin and flavour rotation simultaneously. For Majorana neutrinos the RSFP converts 
left-handed neutrinos of a given flavour into right-handed neutrinos of a different flavour. 
Disregarding flavour neutrino mixing by simplicity the evolution of two Majorana neu- 
trinos {{veLiV^n) in presence of matter and magnetic fields is described by the following 
effective Hamiltonian: 

4(- eL ) = ( V ~J M n Z V- 6 M (4) 



dt \VhrJ V V B ± / V*V? 

where 

Bf(t) = B x (t) ± iB y (t) = \B ± (t)\e^ (5) 

and 

<$>(t) = arctan B y /B x . (6) 

The matter potential is V = y/2Gp(N e — N n ) where N e ,N n are the electron and neutron 
number densities, Gf is the Fermi constant and A = Am 2 /2E. The magnetic field 
strength B enters the evolution Eq.(Q) being multiplied by the neutrino transition moment 
/x. The existing upper limits on the magnetic moment of the electron neutrino include 
the laboratory bound fj, < 2 x 10~ 10 ^£ from reactor experiments as well as stronger (one 
or two orders or magnitude) astrophysical and cosmological limits. Expected values of 
B ~ 1 — 100 kG in the solar convective zone and \x = W~ n /j,B would give an expected 
range for the product fiB ~ 10 -8 — 10 -6 /jgG « 5.6 10 -7 — 10 -15 eV or in the practical 
units which will be used throughout this work \xB « 0.1 — 10 H11B4. 

Assuming that the direction of the transverse magnetic field does not change along the 
neutrino trajectory (absence of twist) the evolution equation can be written as: 

{ ± (v eL \ = / V - A fiB ± \ / u eL \ ^ 



dt \ v^r) V f*B± J \ Vfj,R 

In the case of constant magnetic field and matter density, the corresponding transition 
probability obtained from Eq.(0) is 

PiyeL^V^) = sin 2 (kfiB ± t)/k 2 (8) 

with k 2 = 1 + (A — V) 2 /(fiB±) 2 . In resonance (A ~ V) or for strong fiB± the constant 
k ~ 1 and the probability is a periodic function of B with maximum amplitude ~ 1. 

The problem described by Eq.(|7|) when V is locally varying and B is a constant is 
formally equivalent to a MSW problem with a pseudo mixing angle given by a certain 
function of A and \xB. In the presence of a matter density of exponential profile, as 
approximately happens in the solar case, and a constant transverse field, the precession 



probability can be obtained in a closed analytic form (see [11] for a two generation case 
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|T2|| for a general case) . In this case the evolution matrix for the solution to Eq. (|4|) is given 
by (detailed expressions can be found in fll2|| ): 



[/(Mo) = uUt)U s (t ), U s (t)=exp(iH t)U r (t). (9) 
The Hamiltonian Hq and effective mixing matrix u are 

/ A + \ | / cos 8 - sin 8 \ . . 

H = n n \ n > u = • a a ( 10 ) 

V A_ / cos 8 J 

With A± = A (-1 ± Vl + e 2 ) /2; e = 2fiB_ L /A, tan 2 8 = -A-/A+. The elements of 
the matrix C/ r are simple analytic functions of A, V and the "mixing angle" 8. 

3. The master Equation (||) can be written in terms of the corresponding 2x2 density 
matrix p(t) as: 



= [H reg , P ]+ixB x (t)[a 1 ,p]+ixB y {t){G2, P ]. (11) 



The B x ,B y are the Cartesian transversal components of the chaotic magnetic field. Vac- 
uum mixing terms, matter terms corresponding to the SSM density profile and the regular 
magnetic part Hamiltonian, are all included in H reg . a\^. are the standard Pauli matrices. 

We assume that the components B x ,B y are statistically independent, each of them 
characterized by a <5-correlator. The averaged evolution equation is a simple generalization 
(see H for a complete derivation) of the well known Redfield equation (jH|) for two 
independent sources of noise and reads ( ft 2 = fi 2 {B 2 )L ): 



= [H reg ,{p)]-in 2 {[a l: [a u { P )]] + [a 2 ,[a 2 ,{ P )]}). (12) 



; d(P) 
dt 

Eq.(p^) can be also written as: 

Ap) 



dt 



[H reg , (p)} - 2i£l z (2(p) - axipjax - a 2 (p)a 2 ) . (13) 



In the averaged evolution equations the length Lq appears only through the product tjLq. 
Therefore, we will present our results as a function of the quantity P which is a simple 
function of such a product: 

P = I (1 + exp(- 7 )) , 

7 = AU 2 At = ArjL (pB ) 2 /3. (14) 

One reason for using P is that it is a good approximation for the depolarization that 
the presence of noise induces in the averaged neutrino density matrix. At is the distance 
over which the noise is acting. 

In general Eq.(|l3|) must be solved numerically. For that purpose we have found ad- 
vantageous to reduce such equation to a single integral equation. The procedure, which 
could be of some interest for other applications, is outlined in the appendix. 
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4. We have supposed in our numerical calculation that the noise is effective only in 
a thin layer with thickness At = 0.1 Rq starting at r = 0.65 Rq. Different levels of 
noise are parametrized by the quantity P (Eq.(|i~4|)) which varies between 1 (no noise) and 
1/2 (maximal noise). This equivalent to a range for the r.m.s field ^(B 2 ) ~ — 600 kG 
(assuming the reasonable value of the length Lq = 1000 Km and [i = 10 -11 /i#). 

Fig.(|l|) shows the distortion of the neutrino energy spectrum for different levels of noise. 
We plot the survival probability for an electron neutrino created at the solar interior when 
arriving at the Earth. The neutrino travels in a constant magnetic field of magnitude, 
in this example, [iB = 4.0 H11B4 along all the convective zone. An exponential profile 
as given by [14] is used for the matter density. This plot has been obtained using the 
exact solution (Eqs.@) for the noiseless case and the numerical solution to Eq.(p4|) for 
the non-zero noise cases. Numerical integration of the original differential equations gives 
the same results, at a much more higher CPU cost. 

The squared mass difference used in Fig.(|]) is Am 2 = 5.0 x 10~ 9 eV 2 . This value 
implies that some RSFP occurs at the relevant neutrino energies: For the RSFP to be 
effective in the convective zone it is neccesary that Am 2 /2E < 10~ 7 eV 2 /MeV, or for the 
typical solar energies Am 2 < 10~ 8 eV 2 . It is clear from Fig.(|l|) that the lowest-energy 
pp-neutrinos with characteristic energy E pp 0.3 MeV, which are detected only in the Ga- 
Ge detectors (where they provide more than half of the signal) are hardly affected by the 
magnetic field. Intermediate and higher energy neutrinos are strongly suppressed. When 
noise is applied the distortion tends to disappear and the spectrum becomes practically 
flat over all the energy range considered (for P < 0.7 — 0.8). One important observation 
from the figure is that the influence of the noisy magnetic field is significant for all values 
of Am 2 even for those whose resonance point is not exactly in the region where the chaotic 
field is present. 

In the next figure, Fig.(^), we show the expected signal rates in different neutrino 
experiments. The ratio of the averages of the observed detection rates and the SSM 
expected rates which have been used are shown below (la statistical and systematic com- 
bined errors, [fill ; the SSM predictions are those of [Q): Homestake 0.27 ± 0.027, Su- 
perKamiokande (345 days) 0.375±g;g^, Ga-Ge (SAGE-GALLEX combined) 0.51 ± 0.06. 
For computing the expected signal at the Kamiokande experiment we have taken into 
account the two elastic scattering processes v x e — > v x e with v x = v e ^V ii . We have used the 
expressions for the elastic cross sections appearing in If] with a kinetic energy threshold 
for the observation of the scattered electron of T P 7 MeV. 

We observe that the signal is successively correlated (situation not favored from the 
data) or anti-correlated with the value of the magnetic field. The periodic structure of 
the plots is explained easily if we take into account Eq.(||). The variation for the Ga-Ge 
experiments is always significantly lower than those for Kamiokande or Homestake as it 
could be expected from the previous discussion about the energy spectrum distortion. The 
introduction of noise smoothes the amplitude of the periodic variations 

In Fig.(j3|) we present the single and combined x 2 exclusion plots corresponding to the 

Of course we could also think of time variations induced by a modulated P. 



6 



expected signal rates at the three experiments as a function of the oscillation parameter, 
regular magnetic field and fluctuation amplitude. The Akhmedov (Q) solutions are re- 
produced clearly (Am 2 » 10~ 8 eV 2 , fiB 4 //11-B4 and higher modes). These solution 
regions (shaded areas) do not change with the addition of noise up to a value P = 0.8. 
For the value P = 0.7 we observe that the solution with the lower magnetic field has 
disappeared already. For values P < 0.7 (not shown in the plot) the three allowed regions 
from combined data disappear completely. 

If we observe the outcome of the individual experiments (individual lines) we can con- 
clude that the stronger constraint at high noise comes from the Homestake data. The 
other two experiments show contradictory tendencies which can be explained by the dif- 
ferent weight of the different parts of the neutrino spectrum in each of them. While the 
Kamiokande allowed region tends to shrink with increasing noise, the region allowed by 
the Ga-Ge data in fact enlarges to a maximum extension. At P = 0.7 and low magnetic 
field there is a large common area which is allowed by both Kamiokande and Ga-Ge data. 

Taking into account the three experiments, and within this very simple scenario, par- 
ticle physics solutions to the solar neutrino problem are allowed only if P > 0.8 — 0.7. 
According to the figures, stronger values of the noise are still allowed if we admit solu- 
tions with ever larger magnetic fields. Large magnetic fields in the convective zone are, 
if not completely discarded, at least disfavored by astrophysics. However from the purely 
neutrino experimental side these solutions increase in attractiveness in conjunction with 
large noise because of the (anti)-correlation smoothening observed in Fig.(Q). 

In conclusion we have shown that RSFP solutions to the solar neutrino problem with 
negligible neutrino mixing survive in the presence of random magnetic fields of low or 
moderate level (P > 0.8-0.7, or ^(B 2 ) < 140-200 kG for p = lO" 11 /^ and a fluctuating 
scale Lq = 1000 km). Conversely, by requiring a particle solution to the solar neutrino 
problem, in the proposed scenario this range of values for P can be interpreted as a 
upper limit on the presence small scale fluctuating magnetic fields in the solar convective 
zone. According to ||] the 5-correlation function is a sufficiently good approximation to 
more realistic finite correlators even for relative large correlation length. Moreover, the 
quantity P is connected with the depolarization of the density matrix. It offers some 
advantage with respect to more "physical" parameters as Q 2 or ^(B 2 ): we expect that 
results expressed in terms of P are largely independent of the concrete modelizing of the 
stochastic properties of the fluctuating field (in particular the form of the correlator). 
Additional, potentially stronger, restrictions on the level of noise can be expected when 
considering the experimental bounds in the solar antineutrino flux from Kamiokande and 
the next generation of solar experiments. 
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Appendix 

Our objective in this appendix is to introduce some mathematical manipulations in 
order to obtain an equation adequate for the purpose of numerical calculation. As a 
preliminary step, a term can be eliminated from Eq.(|T3|) if we define a new density (pA(t)) 
as: 

(p) = e- 4 ^-<°V(i)). (15) 

It is convenient in addition to define a vector representation for the density matrix p; we 
define, using the standard Stokes parameters: 

p = Diag ((p u - P22), 25Rpi 2 , 2Qp 12 ) . 

In vectorial form, the equation (Jl^) is: 

Apa) 



1- 



dt 



(fZ^-8tfi 2 Ti)(p A >, (16) 



where and t\ are 3x3 matrices. In particular n, which is vectorial equivalent to 
a\ + o"2, is of the form r\ = i Diag(l, 0, 0). 

Our strategy will be to solve exactly the problem corresponding to the first term of the 
Hamiltonian in Eq.(|i~6|) and to consider the second one as a perturbation. We will be able 
to sum the infinite terms of the perturbative series in terms of a single unknown-function 
to be defined later. This unknown function is the solution of a certain scalar integral 
equation. We have found this procedure, albeit apparently artificious, very advantageous 
in practical computations. 

Let be U°(t,to) the evolution matrix corresponding to the equation: 



= «P (17) 



for the solar matter exponential density profile and for B(t) = Bq, a constant. This 
equation can be solved exactly (see Eqs.([9|-[To|)). Let us consider now the evolution matrix 
U for the full equation Eq.(p^) and its interaction representation counterpart XJj given by: 

U I (t,t ) = U°Ht,t )U(t,t ). (18) 

Our main result, which can be proved easily, is that the elements of Uj are given by the 
following exact formula: 



Uj(t,t ) l3 = 8ij + (-iU 2 ) [ dtiAi(ti)Sj(ti). (19) 
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The vector 5 is the solution of the following integral equation with kernel K: 

S(t) = B(t) + (-iQ 2 ) /* dtiK(t,ti)S(ti). (20) 

Jt 

The vectors A, B are defined as 

Ai(t) = to); Bi(t) = E^(Mo); (21) 

the kernel K is 

K(h,t 2 ) = ^BtitJMh) (22) 
i 

In symbolic form the solution to Eq.(^) can be written as: 

S = (l + itfKy 1 B, (23) 

the solution to Eq.(|l9|) is then, in the same symbolic form, 

Ui = 1 + (-in 2 ) J Ail + itfKy 1 B (24) 

This is the concrete equation which has been used as the basis for numerical compu- 
tations. 
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Figure 1: Distortion of the neutrino energy spectrum for different levels of noise: P = 1.00 
(no noise, solid line); P = 0.95,0.80,0.70,0.50 (consecutive dashed lines). The equivalent 
values for ^/(B 2 ) are respectively 0.0,70,150,200,600 kG (supposing a scale Lq = 1000 
Km, /jl = 10~ n fiB and a noise region width At = 0.1 R®). The corresponding values 
for the ratio r\ are r\ = (B 2 )/Bq = 0,3,14,25,225. It has been used the parameters 
Am 2 = 5.0 10- 9 eV 2 and ^B = 4.0 n u B 4 . 
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Figure 2: The ratio S/Sq (theoretical over observed signal rates) as a function of the 
magnetic field for the three solar neutrino experiments: Combined Ga-Ge (solid line), 
SuperKamiokande (long dashed line) and Homestake (dotted line). Am 2 = 3 x 10~ 8 eV 2 . 
Plot(A), P=l (absence of noise). Plots (B,C,D), P=0.95,0.8,0.7 respectively (see caption 
in Fig.®). 
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Figure 3: Parameter space regions consistent with the observations at the three neu- 
trino experiments. From light to darker shaded areas combined x 2 90,95,99% C.L. inter- 
vals. Lines refer to individual experiments (95% C.L.): Homestake (dotted line), Ga-Ge 
( dot-dashed), Kamiokande (dashed). Plot (A), P=l (absence of noise). Plots (B,C,D), 
P=0.95,0.8,0.7. (See caption in Fig.©). 
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